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An SU (4) model of high-temperature superconductivity and antiferromagnetism has recently been 
proposed. The SO(5) group employed by Zhang is embedded in this 517(4) as a subgroup, suggest- 
ing a connection between our 577(4) model and the Zhang 50(5) model. In order to understand 
the relationship between the the two models, we have used generalized coherent states to analyze 
the nature of the 5*0(5) subgroup. By constructing coherent-state energy surfaces, we demonstrate 
explicitly that the SC/(4) D SO(5) symmetry can be interpreted as a critical dynamical symmetry in- 
terpolating between superconducting and antiferromagnetic phases, and that this critical dynamical 
symmetry has many similarities to critical dynamical symmetries identified previously in other fields 
of physics. More generally, we demonstrate with this example that the mathematical techniques 
associated with generalized coherent states may have powerful applications in condensed matter 
physics because they provide a clear connection between microscopic many-body theories and their 
broken-symmetry approximate solutions. In addition, these methods may be interpreted as defining 
the most general Bogoliubov transformation subject to a Lie group symmetry constraint, thus pro- 
viding a mathematical connection between algebraic formulations and the language of quasiparticle 
theory. Finally, we suggest that the identification of the SO (5) symmetry as a critical dynami- 
cal symmetry implies deep algebraic connections between high-temperature superconductors and 
seemingly unrelated phenomena in other field of physics. 



PACS numbers: 74.20.Mn 

I. INTRODUCTION 

Data for cuprate high-temperature superconductivity 
(SC) suggests d-w&ve singlet pairing in the superconduct- 
ing state and that superconductivity in these systems is 
closely related to the antiferromagnetic (AF) insulator 
properties of the undoped compounds. This proximity of 
superconducting and antiferromagnetic order is unusual 
and suggests that a correct description of cuprate super- 
conductors must permit the SC and AF order to enter 
the theory on a similar footing. 



A. Unified Models 

S.-C. Zhang et al |lj [|, ||] proposed to unify AF and 
SC states by assembling their order parameters into a 
5-dimensional vector and constructing an SO (5) group 
that rotates AF order into d-wave SC order. Recently, 
we introduced an SU (4) model 0] of high-temperature 
SC and AF order having SO (5) as a subgroup. The mo- 
tivation for and methodology of the 517(4) model are 
superficially rather different from those of Zhang and col- 
laborators. However, the appearance of the Zhang SO (5) 
as a subgroup of the S77(4) symmetry implies that the 
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two models must have a strong physical relationship, as 
we have already suggested [El. 

One difficulty in understanding this relationship is that 
the two models are formulated using different approaches 
that employ different languages. In the SU(4) model the 
five operators responsible for AF order (three staggered 
magnetization operator components) and SC order (two 
d-wave pair operators) enter the theory as quantum me- 
chanical operators, on exactly the same footing as all of 
the other ten generators of SU(4) |J. Thus, the SU(4) 
model is a many-body, fully quantum mechanical the- 
ory in which the charge and spin are exactly conserved, 
so there is no spontaneously broken symmetry. This 
is, of course, as it should be: the charge and spin are 
rigorously conserved in a full many-body formulation of 
the problem as they are in nature. In the Zhang SO(5) 
model, the five operators corresponding to the staggered 
magnetization and the d-wave singlet pairing are instead 
treated as order parameters (expectation values of op- 
erators in a broken symmetry state). Therefore, in the 
SO (5) model the antiferromagnetic phase and the super- 
conducting phase are associated with approximate solu- 
tions of the many-body problem that break charge and 
spin symmetry spontaneously. 

The methodology employed in the development of the 
SU(4) model (systematic application of principles of dy- 
namical symmetry) has found broad application in other 
fields such as nuclear structure || 0, molecular physics 
H |{|, or particle physics This implies a deep al- 

gebraic connection between high temperature supercon- 
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ductivity and a variety of phenomena in other fields of 
physics that bear no superficial resemblance to high- 
temperature superconductors. For example, we have 
pointed out previously Q that there is an almost perfect 
mathematical analogy at the algebraic level between the 
AF-SC competition in cuprates and the competition be- 
tween quadrupole deformation of the nuclear mean field 
and nucleon pairing that is a central organizing princi- 
ple of nuclear structure physics. The appearance of the 
50(5) subgroup in the dynamical symmetry of the SU (4) 
model for high-temperature superconductors is very sim- 
ilar to that of the SO (7) subgroup in the SO (8) fermion 
dynamical symmetry model || that describes nuclear 
structure. 



We shall demonstrate here what is, to our knowledge, 
the first example in condensed matter physics of such a 
symmetry. In this case, the critical dynamical symme- 
try will be shown to be based on the SO(5) subgroup 
of SU(4), and it will be shown to interpolate between 
AF and SC order as the hole doping parameter is var- 
ied. Thus, we shall propose that, within the context of 
the more general ST/ (4) model, the SO (5) symmetry em- 
ployed by Zhang serves as a doorway between AF and 
SC order in a manner that can be specified in precise 
terms using the language of Lie algebras and generalized 
coherent state theory, and that is related to critical dy- 
namical symmetries that have been found in other fields 
of physics. 



B. Generalized Coherent States 

There is a well-developed theoretical approach to relat- 
ing a many-body algebraic theory with no broken sym- 
metry to an approximation of that theory that exhibits 
spontaneously broken symmetry: the method of gener- 
alized coherent states pj] |. This method may be viewed 
as the extension of Glauber coherent state theory (which 
is built on an SU(2) Lie algebra) to a more complex 
system having an arbitrary Lie algebra structure. It 
has also been shown to be equivalent to the most gen- 
eral Hatree-Fock-Bogoliubov variational method under 
symmetry constraints, and has been applied extensively 
in various areas of physics and mathematical physics, 
though not to our knowledge in condensed matter. 

The result of such a coherent state analysis is a set 
of energy surfaces that represent an approximation to 
the original theory in which order parameters appear 
as independent variables. In the general case, these en- 
ergy surfaces can exhibit (possibly multiple) minima and 
these minima may appear at non-zero values of the order 
parameters, implying spontaneous symmetry breaking. 
Thus, the generalized coherent state method is a system- 
atic approach to relating a many-body algebraic theory 
to its approximate symmetric and broken symmetry so- 
lutions. 



C. Critical Dynamical Symmetries 

The concept of a critical dynamical symmetry appears 
naturally in applications of generalized coherent state 
techniques to other fields of physics. A critical dynamical 
symmetry is a dynamical symmetry having eigenstates 
that vary smoothly with a parameter (usually particle- 
number related) such that the eigenstates approximate 
one phase of the theory on one end of the parameter 
range and a different phase of the theory at the other 
end of the parameter range, with eigenstates in between 
exhibiting large softness against fluctuations in the order 
parameters describing the two phases. 



D. Symmetry Breaking 

We shall use these results to derive a result that has re- 
ceived considerable attention for the SO(5) model: that 
an exact SO(5) symmetry cannot account for the detailed 
phenomenology of the cuprate superconductors and that 
it is necessary to break SO(5) (explicitly, not sponta- 
neously) in a particular way in order to recover Mott 
insulator normal states at half filling, as is required by 
the data. As we shall show, the embedding of SO(5) as 
a subgroup of ST/(4) implies naturally that SO (5) must 
be broken in this manner in order to produce the correct 
normal states at half filling. 

E. Goals of Paper 

Let us conclude this introduction by enumerating con- 
cisely the primary goals of this paper: 

1. This paper serves to introduce the generalized co- 
herent state method to issues in condensed matter 
physics. 

2. We shall show how to relate the generalized coher- 
ent state to the most general variational quasipar- 
ticle states that can be constructed subject to the 
constraints of SU (4) symmetry. 

3. This paper introduces into the high-temperature 
superconductor discussion in particular and con- 
densed matter in general, the concepts of critical 
dynamical symmetries that have been applied with 
considerable success in other fields of physics. 

4. We shall demonstrate that the coherent state so- 
lution of the SU(4) model identifies SO(5) as a 
critical dynamical symmetry. This critical dynami- 
cal symmetry will be shown to interpolate between 
AF and SC order as the hole doping of the sys- 
tem is varied. This doping dependence is a natural 
consequence of the S£/(4) symmetry, without the 
introduction of a chemical potential ansatz. 
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5. We shall show that the AF and SC phases them- 
selves are more economically described, not by 
50 (5), but by dynamical symmetries built on 
SO (4) and SU(2) subgroups of SU(4), respectively. 
Thus, we shall demonstrate that SU(4) accounts 
for both the origin of the AF and SC order param- 
eters and the 50(5) "rotation" of the superspin 
vector between these phases. 

6. We shall demonstrate that as a fundamental conse- 
quence of the 5C/(4) structure the 5*0(5) subgroup 
must be broken in order to produce Mott insulator 
normal states. Furthermore, we shall demonstrate 
that the required symmetry breaking terms and the 
doping dependence of the solutions occur naturally 
within the 5t/(4) parent algebra and need not be 
introduced by hand as is required in the Zhang 
50(5) model. 

7. We shall show that, because of the nature of the 
critical dynamical symmetry, an AF perturbed 
SU(4:) D 50(5) Hamiltonian is able to approxi- 
mate various symmetry limit solutions depending 
on the doping: the solutions are close to the 50(4) 
limit presenting AF order around half filling, and 
approach the SU(2) limit presenting SC order as 
the hole doping increases. Thus the 51/(4) model 
with a perturbed 50(5) Hamiltonian is able to ac- 
count for the essential features of high T c supercon- 
ductors. 

Our approach will be to introduce the basic features of 
the coherent state technique in sections II and III. In sec- 
tion IV we derive the coherent state energy surfaces for 
the SU (4) model, and in Section V we discuss the SU (4) 
energy surfaces in various dynamical symmetry limits. 
In Section VI, we use SU (4) energy surfaces to examine 
the properties of broken 50(5) symmetry, and then use 
these results in Section VII to argue that with a small 
AF-preferred symmetry breaking, an 5C/(4) D 50(5) 
Hamiltonian may be able to describe high temperature 
superconductivity. Section VIII presents a summary and 
conclusions. 



II. COHERENT STATES AND THE SU(4) 
MATRIX REPRESENTATION 

Gilmore Jl2|, and Perelomov Jl4| (see also earlier 
work by Klauder [^5)) demonstrated that Glauber co- 
herent states for the electromagnetic field could be 
generalized to define coherent states associated with the 
structure of an arbitrary Lie group. In particular, they 
observed that the original Glauber theory for coherent 
photon states may be expressed in terms of an SU{2) 
Lie algebra by examining the commutation properties of 
the second-quantized operators of the theory. Once the 
theory has been formulated in terms of an SU(2) algebra 
generated by combinations of creation and annihilation 



operators, the formalism may be generalized to encom- 
pass a set of such operators closed under any Lie algebra. 

We shall term this extension of the Glauber theory 
to arbitrary Lie algebras the generalized coherent state 
method. Since an extensive review concerning the basic 
approach has been presented in we omit an intro- 
duction to the general technique and proceed directly to 
a specific application of the generalized coherent state 
method to the 5t/(4) algebra. 



A. Algebraic Structure 

A convenient way to analyze these generalized coherent 
states is in terms of their geometry, which is in one-to- 
one correspondence with the coset space. Let us begin 
with the algebraic structure of the SU (4) model . We 
introduce 16 bilinear fermion operators: 



P12 = 5 r ( fc ) c fcT c -fci. Pi2^^9*{k) 



(i) 



'ij ~ 2-^, C k+Q,i C k,j 
k 

J 



where c k i creates a fermion of momentum k and spin 
projection i,j — 1 or 2 = f or |, Q = (tt,tt,tt) is an AF 
ordering vector, f2/2 is the electron-pair degeneracy, and 
following Refs. ^ [| we define 

g(k) — sgn(cosfc :E — cosk y ) 

with the constraints 

g(k + Q) = -g(k), \g(k)\ = l. 

Under commutation the operator set (|l|) closes a t/(4) 
algebra corresponding to the group structure 

D 50(4) x U(l) D SU(2) S x U(l) 
U(4) D SU(4) D 50(5) D SU(2) S x U{1) (2) 
D SU(2) P x SU{2) S D SU(2) S x U(l) 

where we require each subgroup chain to end in the sub- 
group 



SU(2) B x U(l) 



(3) 



representing spin (the SU(2) S factor) and charge (the 
U(l) factor) conservation, because the physical states of 
the system obey these conservation laws. 

In Ref. 0| we discussed the representation structure of 
(||) and showed that the 50(4) subgroup is associated 
with antifcrromagnctism, the SU (2) p subgroup is associ- 
ated with (i-wave superconductivity, and the 50(5) sub- 
group is associated with a transitional symmetry inter- 
polating between the other two. In this paper, we further 
provide the full mathematical justification for interpret- 
ing the 50(5) subgroup as a symmetry interpolating dy- 
namically between SC and AF phases. 
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B. A Convenient Basis for Generators 



that the following mapping preserves the algebra of i 



It is convenient to take as the generators of U(4) 
[/(l) c d x SU(4) the new combinations 

0+ = Qn + Q22 = X^ c fe+QT Cfe T + C k+Ql c kl) 



S = 
Q = 

#t - 



$12 + $21 ■ S\2 — S21 Sn — S22 

= , -I = , = 



Q12 + Q21 . Q12 — Q21 Qn — Q22 
, —1 , 



(4) 



• 9ii Q22 I11 + Q22 



■912+921 



I 2 ' 2 ' 2 I 

7? = (7?t)t D^=p\ 2 D= Pl2 M= 1(511 +5a2) 

where Q+ generates the J7(l) c d factor and is associated 
with charge density waves (do not confuse this U(l) fac- 
tor with the one appearing in Eq. (||) that is associated 
with charge conservation), 5 is the spin operator, Q is 
the staggered magnetization, the operators 7?^ and if are 
those of Rcf . pj , the operators and D are associated 
with d-wave pairs, and 

2M = n - 9. 

is the charge operator. Because of the direct product 
structure £/ (4) — > f/(l) c d x 517(4), we can without loss 
of generality analyze the t7(4) structure in terms of its 
subgroup 5t/(4), with the /7(l) cc j factor considered sep- 
arately. Hence all subsequent discussion will deal with 
the 5£/(4) subgroup of U{A). 

To facilitate comparison with the 5*0(5) symmetry, the 
50(6) ~ SU(A) generators may be expressed as 
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(5) 



where we define 

D± = \(D±D^) ^± = 1(^+4) (6) 
with 7 a & = — and with commutation relations 

[L a b, L c d] = i(5 ac Lbd — SadLbc — $bcL a d + <5bdL a c)- (7) 
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where 


0^ = 


— Qii + 2"? 





(8) 



the a x , <7 V , and a z are Pauli matrices in the standard 
representation, 

a± = \{cr x ± ifj/)' 

and J is a unit matrix. Likewise, it is easily verified that 
in terms of this representation, 
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C. Faithful Matrix Representation 



D. Collective Subspace 



The coherent state method requires a faithful matrix 
representation of 5J7(4). Explicit multiplication verifies 



We take as a Hilbert space 
|S) = \n x n y n z n d ) = (4)^(^(4)^(0^0) 
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which is a collective subspace associated with 5*0(6) ir- 
reps of the form 



0,0). 



In this notation we use the well-known isomorphism of 
SU(4) and 5*0(6) to label the irreducible representations 
with the standard 50(6) quantum numbers (ci, 02, 03) 
IPTf . Physically, this irrep represents a "maximally 
stretched" state in the representation space that is in 
turn associated with maximal collectivity; as such, it is 
the obvious candidate for a collective subspace describing 
the lowest states of the system. 

One sees immediately that the expectation value of Q+ 
is zero for any state in this collective SU(4) space: the 
matrix representation of 0+ is a un it matrix, and thus 
Q + commutes with all the SU(4) generators, leading to 

(5|Q+|5)=0. 

In the symmetry limit, this implies that charge density 
wave excitation are excluded from the 5C/(4) model re- 
stricted to this subspace |lq|. 



E. SU(4) Casimir Operator 

The Casimir operator of 5Z7(4) - SO (6) 
C su{ a) = 7? f -7? + £>t D + S • 5 + Q ■ Q + M(M - 4) (10) 

is an invariant and its expectation value in this collective 
subspace 



(11) 



is a constant. 



III. SYMMETRY-CONSTRAINED 
BOGOLIUBOV TRANSFORMATION 

Utilizing the methods of Ref. , the coset space is 
SU{A)/SO{A) x 1/(1), 

where the 50(4) subgroup is generated by Q and 5, and 
[7(1) is generated by the charge operator M. The coher- 
ent state may be written as 

|V>=T|0*>. (12) 

The operator T is defined by 

T = exp(r?ooPi 2 + VwqU ~ h. c), (13) 

where |0*) is the physical vacuum (the ground state of the 
system), the real parameters 7700 and 77x0 & re symmetry- 
constrained variational parameters, and h. c. means the 
hermitian conjugate. Since the variational parameters 



weight the elementary excitation operators p\ 2 and q\ 2 
in Eq. (|l3|), they represent collective state parameters 
for a subspace truncated under the 5t/(4) symmetry. 
The most general coherent state corresponds to a 4- 
dimensional, complex, compact manifold parameterized 
by 8 real variables. The reduction of the coherent state 



parameters to only two in Eq. (13) follows from requiring 
time reversal symmetry and assuming the conservation of 
spin projection S z for the wavefunction. 

It is often simpler to view the coherent states as 
Hatree-Fock-Bogoliubov (HFB) variational states con- 
strained by the dynamical symmetry. The symmetry- 
constrained HFB coherent state method is discussed in 
Refs. @ |0|, §3, It may be viewed as a type of 
mean-field approximation to the underlying many-body 
problem that is particularly useful in the present context 
because it leads to easily visualized energy surfaces. This 
identification provides a natural connection to sponta- 
neously broken symmetries and effective Lagrangian field 
theories on the one hand, and to quasiparticle language 
on the other. 

From the coset representative expressed in the 4- 
dimensional matrix representation (|J), the transforma- 
tion operator T defined in Eq. ( p"2| ) may be written as 



T 



Y 1 X 
-Xt Y 2 



X 



a + (3 
-{a- (3) 



where Yi and Y 2 are determined by the requirement 
that T be unitary, and a and are variational param- 
eters related to 7700 and 7710 in Eq. Jlg| ) (see Ref. 
Introducing 



v + = a + v_ = a — (3, 
the requirement of unitarity gives 



(14) 



X 





— v_ 







with the constraint that 
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u+ 

u_ 



vi = l 



u_ 
u + 



(15) 



A. Quasi-fermion Transformation 



The existence of the 4-dimensional matrix represen- 
tation of the SU(4) algebra implies the existence of a 
representation in which the single-particle basis can be 
written in the form 



where a is a state conjugate to a. In this representation, 
any generator O may be written as 



Qll Q12 

021 22 



o 



E 



(12) c f .c f ■ 



[0 (11) c t c 

ol 21) . 



(22) c-c f 
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and Eq. (12) is seen to be a Bogoliubov type transforma- 
tion, but one that is constrained to preserve the SU(4) 
symmetry. 

Through the operator T, the physical vacuum state 
|0*) (the ground state of the system) is transformed to 
a quasiparticle vacuum state with the parameters a 
and (3 (or v±) determined by minimizing the energy of 
the system. Likewise, the basic fermion operators 

{ C !tT' C <t|' Cff T' c &i\ 

are transformed to quasifermion operators 



through 



T 



\4J 



(16) 



B. One-Body and Two-Body Operators 



Using the transformation (|16|), one can express any 
one-body operator in the quasiparticle space as [E3| 



TOT- 1 



£,(11) 0(12) 
0(21) 0(22) 



o = E o « w + EK )a Ui-o},'«ii0« 



i( 22 )„t 



+ (1 V a f 

T ^i.j a cri a 3j 



where the 's are fixed by the transformation prop- 



O 



21) \ 



erties of the operator O: 



O 



{ixv) 



_ <j(lim) Q(mn) — l^ni/Jj^ 



(17) 



and T^ m ) and C>( mri ) are two-dimensional submatrixes 
of T and O, respectively. 

Because the quasiparticle annihilation operator acting 
on the quasiparticle vacuum \ip) is zero, the expectation 
values for one-body operators O are given by 

(6) = MOM =5>i 22) =£Tr(0( 22 )) (18) 



and for two-body operators Oa Ob, 

(d A o B ) = (tP\6 a o b \tP) 

= ^Tr^f^Tr^) 



(19) 



C. Expectation Value of Generators 

Utilizing equations (|9|) and ()l^)-(|l9|), and noting that 
the summation provides a factor of fi/2 since the 
matrix elements of Eq. (|l^) do not depend on a, one 
obtains the expectation values for all the generators and 
their scalar products in the coherent state representation: 



(Z>t) = 


(D) = 


if2(u_v_ + u^ 


■v+), 


(20) 


(4) = 


(T«> = 


if2(u_v_ — u. 


f-v+), 


(21) 


Q = 


(G*> = 


= ±OK-v 2 _) 




(22) 


(n) = 


n = 


^v^ + v 2 .), 




(23) 


(fx) = 


K) = 


(S) = (Q*) = 


(e«) = o, 


(24) 


(DtD) = 


If7 2 (u 


f v + + u_v_) 2 


+ ifi(v| ^ 


-v 4 _), (25) 


7T'-7r) = 


I^ 2 (u 


fV+ - u_v_) 2 







+ ±fi(<+v 4 _+4v 2 V), 

;q- q) = ^ 2 (v 2 h -v 2 _) 2 + ^[(u + v + ) 2 

+ (u_v_) 2 + (u + v_) 2 + (u_v + ) 2 ; 
[~S-^) = ifi[(u + v_) 2 + (u_v + ) 2 ], 
(M 2 ) 



(26) 

(27) 
(28) 

i(n-f7) 2 + if7[(u + v + ) 2 + (u_v_) 2 ]. (29) 

Using the above results, one can also verify Eq. (|ll|) ex- 
plicitly for the expectation value of the Casimir operator 

Cs«(4) ■ 



D. Order Parameters 

By virtue of the unitarity condition ([l5|), there are 
only two independent variational parameters in the above 
equations. They may be chosen as either v + and v_ , or as 
a and j3, using fll4|). However, from Eq. (|2^) the squares 
of v± (or of a and (3) are constrained by the equation of 
a circle since 



(n) = n{v 2 + + v 2 .) = 2n(a 2 + (3 2 ). 



(30) 



Thus, for a fixed particle number n we may evaluate ma- 
trix elements with only a single variational parameter, 
say (3, which may in turn be related to standard order 
parameters by comparing matrix elements. For exam- 
ple, the z component of the staggered magnetization is 
related to (3 and v± by 

Q = (Q z ) = ±r!(v+-v 2 _) 

= 2ft/3(ra/(2fi) -/3 2 ) 1/2 . (31) 

These measures of antiferromagnetic order are in turn re- 
lated to the superconducting order parameter a through 
Eq. (§(]). From Eqs. @ and Q, the ranges of (3 and 
a are 



< (3 < vV 40 < ol < y/n/2tt. 



Using Eqs. (|20|)-(|3l|) one can then evaluate the energy 
surface as a function of Q or (3 or a, and study the ground 
state properties of the SU (4) model. 
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IV. COHERENT STATE ENERGY SURFACES 



The most general Hamiltonian for the S£/(4) model 



H = en 



xQ-Q + gS-S, 



(32) 



where e, v, Go, G\, x> and g are parameters defining the 
strengths of single-particle and interaction terms. Since 
C su (4) is an ST/ (4) invariant, if we assume for the ground 
state spin that (S) — and that the number of particles 
n is a good quantum number, Eqs. (flO|)-(P"l|) imply that 
the Hamiltonian (J3^) may be parameterized without loss 
of generality as 

H = H - G [ (1 -p)D J< D +pQ ■ Q ], (33) 
where p lies in the interval to 1 and 



G 



(0) 
off 



(l-p)G = G -G 1 



Xcff = pG = X ~ Gi 
Ho = en — vn 2 — G\ 



(34) 



l-x 2 



and Xeff ar e the effective strengths of pairing and 
the <2 ■ Q interactions, respectively. 
From (^0|)-(^l|), one can show that 



v 2 = — ± — 
V± 2VL Q ' 

Eqs. (p0|) and (Blh can be written as 



(35) 



A = (Z?t) = (_D) = Vz?Fd = A+ + A_, (36) 
n = (ttJ) = (tt 2 > = Vyft - tt = A+ - A_ (37) 



where 



A+ = i^A 



T 2 



and x is the effective hole concentration 

n 



x = l- 



0' 



(38) 



(39) 



(By "effective" we mean that X is a ratio of the hole-pair 
number (£1 — n)/2 to the pair degeneracy 0/2, rather 
than the ratio of hole number to the total number of 
lattice sites.) It can be estimated that to avoid hole- 
pair collapse, O is required to be roughly one-third of the 
total lattice sites, and in turn the true hole concentration 
is one-third of x Q . 

The quantities A and n present the spin-singlet and 
spin-triplet pairing correlations. The former is propor- 
tional to the superconducting pairing gap; the latter can 



be regarded as a measure of the SO(5) correlation since 
the 5*0(5) Casimir operator is 



C, 



sob 



TT - TT + S ■ S + M{M -3). 



(40) 



By utilizing Eqs. (P5Q-(]27|), a general expression for the 
energy surface of the SU(A) Hamiltonian as a function 
of the antiferromagnetic order parameter Q may be ob- 
tained. In the Q — > oo limit, the energy surface is defined 
by 



E{Q) = (H) -H = -Go [(1 - P )A 2 + PQ 2 



(41) 



Converting Q into the alternative order parameter [3 al- 
lows us to express the energy surface as a function of /3 
and n, 

G n 2 



2 (n_ 




i _^y_w^ 




which may also be expressed in terms of the supercon- 
ducting order parameter a using (J3C 



COHERENT STATES AND SU(4) 
SUBGROUPS 



Assuming Go > (suggested by phenomenology), 
p = 1/2 in Eq. j33| ) corresponds to SO(5) symmetry, 
while the extreme values and 1 correspond to SU(2) 
and 50(4) symmetries, respectively (see Ref. Q). Other 
values of p respect 5C/(4) symmetry but break the 50(5), 
50(4), and SU(2) subgroups. In Fie^l we illustrate the 
ground-state energy E(0) of Eq. ( |42j ) as a function of 
the order parameter (3 for different electron occupation 
fractions n/VL with p = 0,| and 1. 



A. SU(2) Limit 

For p = [SU(2) limit; see Fig. la], the minimum 
energy occurs at j3 = (equivalently, Q = 0) for all 
values of n. Thus, A reaches its maximum value of 



(43) 



A — ifK/1 - r 2 
indicating superconducting order. 

B. SO(4) Limit 



For p — 1 [50(4) limit; see Fig lc], the opposite sit- 
uation occurs: f3 = is an unstable point and an in- 
finitesimal fluctuation will drive the system to the energy 
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FIG. 1: Coherent state energy surfaces. The energy units are G^fi 2 /4 for figures (a) and (b), and Xcff^ 2 /4 for (c). Ho is 
taken as the energy zero point. Numbers on curves are the effective lattice occupation fractions, with n/fi = 1 corresponding 
to half filling and < n/tt < 1 to finite hole doping. 5*0(5) symmetry corresponds to p = 1/2 and the allowed range of j3 
is [—\\J n/fi, \ y/n/Q], which depends on n. The order parameter f3 is related to the order parameter Q = (Q z ) (staggered 
magnetization) and the electron number n through Eq. (plj). 



minima at finite, 



D. SO(5) As a Critical Dynamical Symmetry 



[3 = ±±^/n. 



(44) 



Thus, \Q\ reaches its maximum value of n/2, indicating 
the presence of AF order. 



C. SO(5) Limit 

From Fig. lb, the 50(5) dynamical symmetry (with 
p — i) is seen to have extremely interesting behavior: the 
minimum energy occurs at (3 = Q for all values of n, as 
in the SU{2) case, but there are large-amplitude fluctu- 
ations in AF and SC order. In particular, when n is near 
Q, (half filling), the system has an energy surface almost 
flat for broad ranges of (3 (or Q or a). This suggests a 
phase having much of the character of a "spin glass" for a 
range of doping fractions (specifically, a phase very soft 
against fluctuations in the order parameters). It could 
also lead to inhomogeneous structures such as stripes 
if there is a periodic spatial modulation of the system, 
since the soft nature of the energy surface implies that 
relatively small perturbations can shift an SO (5) system 
between antifcrromagnctic and superconducting behav- 
ior. As n/ft decreases, the fluctuations become smaller 
and the energy surface tends more and more to the SU (2) 
(superconducting) limit. 



Dynamical symmetries that, within the dynamical 
symmetry itself, exhibit a transition between qualita- 
tively different energy surfaces as a parameter (usually 
related to particle number) is varied have been termed 
critical dynamical symmetries The 5*0(5) dynam- 

ical symmetry, within the context of its SU(4) parent 
symmetry, exhibits such transitional properties. At half 
filling the energy surface is completely flat under varia- 
tions of the antiferromagnetic order parameter (3 (see the 
n = 1.0 curve of Fig. lb), implying large fluctuations in 
the order parameters. But as hole doping is increased the 
50(5) energy surface changes smoothly into one localized 
around (3 — (see the n = 0.1 curve of Fig. lb). Thus, 
50(5) is an example of a critical dynamical symmetry. 
Such symmetries are well known in nuclear structure 



physics 1 11, |19|, |20(]. The 50(5) critical dynamical sym- 
metry discussed here in a condensed matter context has 
many formal similarities with the 50(8) D 50(7) critical 
dynamical symmetry of the (nuclear) Fermion Dynami- 
cal Symmetry Model H. The condition for realization 
of the 50(5) critical dynamical symmetry is that the 
strength of Q ■ Q equals that of D'D in the Hamilto- 
nian; This is similar to the SO (7) nuclear critical dynam- 
ical symmetry, which is realized when there is an over- 
all 50(8) symmetry and the monopole pairing and the 
quadrupole interaction terms are of equal strength in the 
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nuclear Hamiltonian In the SO(7) case of nuclear 

physics, the order parameter analogous to (3 presents nu- 
clear deformation: nuclei around midshell (half filling of 
a shell by nucleons) are soft against shape fluctuations 
and transform into a spherical shape (the 50(5) x SU(2) 
dynamical symmetry limit of the 50 (8) symmetry) as 
the number of nucleons increases. 



variations of the AF, SC and S0(5) correlations (Q, A 
and II) with the hole doping x are shown in Fig. 2c, while 
the variations of the contributions of each term in the 
Hamiltonian (the pairing D'D and the AF interaction 
Q ■ Q) to the total energy are shown in Fig. 2d. In Fig. 
2, there is an important quantity, the critical doping x c , 
which can be expressed analytically as |e3 



VI. SO(5) SYMMETRY BREAKING 

Under exact SO (5) symmetry (p = i), the AF and SC 
states are degenerate at half filling. There is no barrier 
between AF and SC states, and one can fluctuate into the 
other at zero cost in energy (see the n/Q = 1 curve of 
Fig. lb) . This situation is inconsistent with Mott insulat- 
ing behavior at half-filling. The Zhang 50(5) model has 
been challenged because under exact symmetry it does 
not fully respect the phcnomcnological requirements of 
"Mottness" . As Zhang H has recognized, for antiferro- 
magnetic insulator properties to exist at half filling, it is 
necessary to break 50(5) symmetry. Such breaking of 
the 50(5) subgroup symmetry is implicit in the £{7(4) 
model, occurring naturally in the 5/7(4) model if p > 1/2 
in the Hamiltonian (^3|). Furthermore, the 5/7(4) sym- 
metry leads to the following constraint 



x r . = \ 1 - 



G 



(0) 
eg 



(D^D + Q-Q + if-n) 



1 



(45) 



This ensures a doping dependence in the solutions, which 
is necessary for describing the transition from AF to SC 
in the cuprates. 

Thus, the coherent state analysis indicates that the 
phenomenologically required 50(5) symmetry breaking 
and the doping dependence in the solutions occur natu- 
rally in the SU (4) model. They need not be introduced 
empirically as proposed in the original Zhang 50(5) 
model 0. Recently, a projected 50(5) model j25| has 
been introduced. Its essence is a patch to the origi- 
nal 50(5) model that implements the Gutzwiller pro- 
jection in order to satisfy the large-/7 Hubbard (non- 
double-occupancy or the Mott-insulator) constraint. In 
our 5/7(4) model, there is no need to introduce such 
a projection artificially because the 5C7(4) symmetry 
constraint already implies a constraint of non-double- 
occupancy with charge density localized on sites of the 
underlying lattice. We shall demonstrate this explicitly 
and give a detailed discussion of the consequences in a 
subsequent paper. 

To see in more detail how in the SU (4) model a broken 
50(5) symmetry can interpolate between AF and SC 
states as particle number varies, let us perturb slightly 
away from the 50(5) limit of p = 1/2 in Eq. @. In Fig. 
2a, 5/7(4) coherent state results for p = 0.52 are shown. 
We denote the value of f3 minimizing (H) as /3q. The 
corresponding variation of the AF order parameter Q = 
(Qz) with n, and its comparison with the variation in 
various symmetry limits are summarized in Fig. 2b. The 



= a/1- 



1-p 



XeS 



For p = 0.52, we have x c = 0.277. 



(46) 



A. Antiferromagnetic Order 

One sees from Fig. 2 that if n is near (half filling) , 
(3q ~ ±0.5; this corresponds to AF order, since the stag- 
gered magnetization reaches its maximum, Q — f2/2, and 
there are no pairing or 50(5) correlations (A — H = 0). 
With the onset of hole doping, n/Q decreases (x in- 
creases). The AF correlation Q quickly diminishes and 
the pairing and 50(5) correlations, A and n, increase. 



B. Underdoped SO(5) Fluctuations 

Before Q vanishes at the critical doping x c — 0.277 
(n/fl — 0.723), the system has an energy surface almost 
flat for broad ranges of (3, implying the presence of large- 
amplitude fluctuations in AF order (and equivalently in 
SC order). Meanwhile, the 50(5) correlation n increases 
and reaches its maximum at the doping where the pairing 
and AF correlations become equal to each other (see Figs. 
2c and 2d). This is the underdoped SC region. The 
coexistence of these three correlations competing with 
each other is consistent with the complexity and variety 
of experimental phenomena in this region. 



C. Superconducting Order 

For small values of n (x > x c ), the stable point is 
/?o = 0. This corresponds to SC order, since both the AF 
and 50(5) correlations vanish (Q = II = 0), and only the 
pairing correlation remains (A > 0). The critical doping 
x c is the optimal doping point since A is maximum at x c 
and decreases as hole-doping increases. Thus the doping 
range x > x c may be considered to be the overdoped 
region. The critical doping x c depends on the ratio of the 
pairing and the Q - Q strengths (see Eq. ((46|)). The larger 
the pairing strength G^ relative to the Q ■ Q strength 
Xeff in Eq. (Eg), the smaller the critical doping value x c . 
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FIG. 2: (a) As for Fig. 1, but for slightly perturbed SO(5). The dotted line indicates the location of the ground state as 
n varies and the corresponding f3 value is denoted as /3q. (b) Variation of the AF order parameter with effective occupation 
number for different values of p. (c) Variation of the AF, SC and SO(5) order parameters Q, A, and II as functions of the 
effective hole concentration x. (d) Variation of the ratio of pairing and Q ■ Q interactions to the total energy of the system as 
functions of the effective hole concentration x. 



VII. HAMILTONIAN FOR HIGH T c 
SUPERCONDUCTIVITY 

From the above discussion one can see that, with a per- 
turbed SO(5) symmetry, the system can undergo phase 
transitions from the AF order at half filling to the SC 
order at smaller filling as particle number varies. This 
picture is at least qualitatively consistent with the obser- 
vations. The 50(5) symmetry breaking in the Hamilto- 
nian (p larger than 1/2) is crucial. Only when 5*0(5) is 
broken does the energy surface interpolate between AF 
and SC order as doping is varied (compare the surfaces 
for p = 1/2 and p = 0.52 in Figs. 1 and 2). We thus 
conclude that high temperature superconductivity may 
be described by a Hamiltonian that conserves ST/ (4) but 
breaks (explicitly) 50(5) symmetry in a direction favor- 
ing AF order over SC order. 

The deviation from the 50(5) symmetry need not be 
large. Experimentally, it is known that the optimal dop- 
ing P c is around 0.16, suggesting that x c ~ 0.48 (note 
that x c ~ 3P C ). This leads to p = 0.56 according to Eq. 



(|46|), which is formally quite close to the SO (5) symmetry 
limit. 

However, it should be stressed that because of the crit- 
ical nature of the SO(5) symmetry, a slightly perturbed 
SO (5) Hamiltonian may have solutions that are close to 
the other symmetry limits of the St/ (4) model for partic- 
ular electron occupation ratios, even though the Hamil- 
tonian itself is not formally in any of the dynamical sym- 
metry limits. For example, one can see from Figs. 2c 
and 2d that the Hamiltonian near half filling actually 
behaves like an SO(4) Hamiltonian since it effectively 
contains only the Q ■ Q correlation term (Q ■ Q is the 
primary component of the SO (4) Casimir). Likewise the 
perturbed SO(5) Hamiltonian approximates the St/(2) 
Hamiltonian containing only the D> D term when x > x c 
(see Figs 2c and 2d) . Only in the intermediate doping 
range (0 < x < x c ), where both D^D and Q ■ Q corre- 
lations have significant contributions, does the p = 0.52 
Hamiltonian behave like an approximate SO(5) Hamilto- 
nian, as it should formally. In particular, near the region 
where the D'D and Q - Q terms have equivalent contribu- 
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tions and the 50(5) correlations (the operator II) reach 
their maximum, the p — 0.52 solution lies very close to 
the 50(5) symmetry limit, as one would expect. 

The present analysis implies that the underdoped 
regime is naturally associated with the 5(7(4) D 50(5) 
dynamical symmetry interpolating between antiferro- 
magnetic and superconducting order. Likewise, opti- 
mally doped and overdoped superconductors are natu- 
rally associated with the SU (4) D SU (2) dynamical sym- 
metry and AF insulators near half filling are associated 
with the 5(7(4) D 50(4) dynamical symmetry (or small 
perturbations around these symmetries). 

As we shall discuss in a separate publication ]2q ], the 
appearance of pseudogap behavior |27] can be described, 
and much of the quantitative phase diagram in cuprates 
can be reproduced rather well p4| , if a fixed Hamilto- 
nian with slightly broken 50(5) symmetry but preserving 
5(7(4) overall symmetry is adopted. This again supports 
our interpretation of 50(5) symmetry as an critical dy- 
namical symmetry. The small 50(5) symmetry breaking 
distorts the completely flat energy surface at half filling 
to stabilize AF character in the system; the critical na- 
ture of the 50(5) dynamical symmetry (that the system 
interpolates between two phases) remains. 

VIII. SUMMARY AND CONCLUSIONS 

The present paper serves to introduce into the con- 
densed matter literature the technology of generalized 
coherent states. As we have shown in the example dis- 
cussed here, these methods provide a systematic way to 
relate a many-body theory to its approximate broken 
symmetry solutions. This approach may be viewed as a 
standardized technology for constructing energy surfaces 
for many-body theories defined in terms of the algebra of 
their second-quantized operators, or equivalently as the 
most general Bogoliubov transformation permitted, sub- 
ject to a symmetry constraint on the Hamiltonian of a 
system. 

To illustrate the power of these techniques, we have 
used generalized coherent states to understand the rela- 
tionship between the SU (4) model of superconductivity 
@ and the Zhang 50(5) model §. The use of 5*7(4) co- 



herent states to analyze the energy surface of its 50(5) 
subgroup permits us to interpret the 50(5) as a criti- 
cal dynamical symmetry that interpolates between AF 
and d-wave SC order as doping is varied, and suggests 
similarities with analogous critical dynamical symmetries 
well known from nuclear structure physics. This permits 
the 50(5) symmetry to be understood dynamically as 
a critical phase that, for a range of doping, has an en- 
ergy surface extremely soft against AF fluctuations and 
therefore having much of the character of a spin glass 
(or possible stripe phases in the presence of a spatially 
modulated perturbation) . 

Thus, the coherent state analysis suggests that 50(5) 
is the appropriate symmetry of the underdoped regime, 
but that the AF phase at half filling and the optimal and 
overdoped SC phases are described by two other 5(7(4) 
subgroups: 50(4) and SU(2), respectively. The coher- 
ent state analysis also shows clearly that the requirement 
of small deviations from 50(5) symmetry and the neces- 
sary doping dependence of the solutions that are inserted 
in the Zhang model occur naturally when 50(5) is a sub- 
group of 5(7(4). 

In addition, we note that the results obtained here may 
have some broader implications. Although the present 
application is specifically to the high-temperature super- 
conductor problem, we may anticipate that these math- 
ematical techniques could find use for any application in 
condensed matter physics where it is important to un- 
derstand the relationship between an exact many-body 
theory and the order parameter(s) characterizing its ap- 
proximate broken symmetry solutions. Clearly there are 
many such possibilities. 

Finally, the concept of a critical dynamical symmetry 
that we have introduced here in a condensed matter con- 
text is one that has already found important application 
in other areas of physics. This implies that there may 
be deep algebraic analogies between various condensed 
matter systems and superficially different systems ap- 
pearing in other fields of many-body physics. We have 
suggested one such analogy here between the physics of 
high-temperature superconductors and the physics of col- 
lective states in heavy atomic nuclei. 
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ural Science Foundation of China. 
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